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ABSTRACT 

Even though the existence of intermediate-mass black holes (IMBHs, black holes 
with masses ranging between 1O 2 " 4 M ) has not yet been corroborated observationally, 
these objects are of high interest for astrophysics. Our understanding of the formation 
and evolution of supermassive black holes (SMBHs), as well as galaxy evolution mod- 
eling and cosmography would dramatically change if an IMBH were to be observed. 
From a point of view of traditional photon-based astronomy, which relies on the mon- 
itoring of innermost stellar kinematics, the direct detection of an IMBH seems to be 
rather far in the future. However, the prospect of the detection and characterization of 
an IMBH has good chances in lower-frequency gravitational- wave (GW) astrophysics 
using ground-based detectors such as LIGO, Virgo and the future Einstein Telescope 
(ET). We present an analysis of the signal of a system of a binary of IMBHs (BBH from 
now onwards) based on a waveform model obtained with numerical relativity simula- 
tions coupled with post-Newtonian calculations at the highest available order. IMBH 
binaries with total masses between 200 — 20000 M & would produce significant signal- 
to-noise ratios (SNRs) in advanced LIGO and Virgo and the ET. We have computed the 
expected event rate of IMBH binary coalescences for different configurations of the 
binary, finding interesting values that depend on the spin of the IMBHs. The prospects 
for IMBH detection and characterization with ground-based GW observatories would 
not only provide us with a robust test of general relativity, but would also corroborate 
the existence of these systems. Such detections should allow astrophysicists to probe 
the stellar environments of IMBHs and their formation processes. 
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1. Motivation 

By following the stellar dynamics at the center of our Galaxy, we have now the most well- 
established evidence for the existence of a SMBH. The close examination of the Keplerian orbits 
of the so-called S-stars (also called SO-stars, where the letter "S" stands simply for source) has 
revealed the nature of the central dark object located at the Galactic Center. By following S2 
(S02), the mass of SgrA* was estimated to be about 3.7 x 1O 6 M within a volume with radius no 
larger than 6.25 light-hours ( [Schodel et al.||2003|[Ghez et al.||2003| ). More recent data based on 16 



years of observations set the mass of the central SMBH to ~ 4 x 1O 6 M ( Eisenhauer et al.||2005 
Ghez et al.|2005|[2008||Gillessen et al.|2009l ). 



Massive black holes in a lower range of masses may exist in smaller stellar systems such as 
globular clusters. These are called intermediate-mass black holes (IMBHs) because their masses 
range between M ~ 10 2 and M ~ 1O 4 M , assuming that they follow the observed correlations 
between SMBHs and their host stellar environments. Nevertheless, the existence of IMBHs has 



never been confirmed, though we have some evidence that could favor their existence (see Miller 
& Colbert|2004l|Miller|2009| and references therein). 



If we wanted to apply the same technique to detect IMBHs in globular clusters as we use for 
SMBHs in galactic centers, we would need ultra-precise astronomy, since the sphere of influence 
of an IMBH is ~ few arc seconds. For instance, for a 10 4 A/ IMBH, the influence radius is of ~ 5" 
assuming a central velocity dispersion of a = 20kms -1 and a distance of ~ 5 kpc. The number 
of stars enclosed in that volume is only a few. Currently, with adaptive optics, one can aspire - 
optimistically- to have a couple of measurements of velocities if the target is about ~ 5 kpc away 
on a time scale of about 10 yrs. The measures depend on a number of factors, such as the required 
availability of a bright reference star in order to have a good astrometric reference system. Also, 
the sensitivity limits correspond to a K-band magnitude of ~ 15, (B- MS stars at 8 kpc, like e.g. 
S2 in our Galactic Center). 

This means that, in order to detect an IMBH or, at least, a massive dark object in a glob- 
ular cluster center by following the stellar dynamics around it, one has to have recourse to the 
Very Large Telescope interferometer and to one of the next-generation instruments, the VSI or 



GRAVITY ( |Gillessen et al.||2006[ |Eisenhauer et al.||2008[ ). In this case we can hope to improve 



the astrometric accuracy by a factor of ~ 10. Only in that scenario would we be in the position of 
following closely the kinematics around a potential IMBH so as to determine its mass. 

GW astronomy could contribute to IMBH detection. In the past years, the field has reached 
a milestone with the construction of an international network of GW interferometers that have 
achieved or are close to achieving their design sensitivity. Moreover, the first-generation ground- 
based detectors LIGO and Virgo will undergo major technical upgrades in the next five years that 
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will increase the volume of the observable universe by a factor of a thousand^] 

The data that will be taken by the advanced interferometers are expected to transform the 
field from GW detection to GW astrophysics. The availability of accurate waveform models for 
the full BBH coalescence in order to construct templates for match-filtering is crucial in the GW 
searches for compact binaries. The construction of this kind of templates has recently been made 
possible thanks to the combination of post-Newtonian calculations of the BBH inspiral and nu- 
merical relativity simulations of the merger and ringdown. Two approaches to this problem are the 
effective-one-body techniques (Buonanno & Damour 1999; [Buonanno et al. 2009| ) and the phe- 
nomenological matching of PN and NR results ( [Ajith et aL"1|2008[ |2009| |Santamaria et al.||20T0| ). 
In this paper we use the results of the latter. 

The structure of this paper is as follows: In section [2] we expand the astrophysical context 
to this problem and give a description of the different efforts made to address the evolution of 
a BBH in a stellar cluster from its birth to its final coalescence. In section |3] we introduce the 
techniques used in the data analysis of the waveform modeling of BBH coalescences, present our 
hybrid waveform model, and use the model to compute and discuss expected signal-to-noise ratios 
in present and future GW detectors. The use of the new waveform model allows us to give an 
improved estimate for the number of events one can expect for several physical configurations in 
Advanced LIGO and the Einstein Telescope, which we present in section |4} We conclude with a 
summary of our results and future prospects of our work in section [5] 



2. Life of a massive binary 

The aim of this section is not to give a detailed explanation of the processes of formation of 
IMBHs and BBHs, but a general introduction of the two different scenarios that play a role in the 
formation of BBHs. 



2.1. Birth 

Up to now, the IMBH formation process which has drawn the most attention is that of a young 
cluster in which the most massive stars sink down to the center due to mass segregation. There, 
a high-density stellar region builds and stars start to physically collide. One of them gains more 
and more mass and forms a runaway star whose mass is much larger than that of any other star in 
the system, a very massive star (VMS). Later, that runaway star may collapse and form an IMBH 



http://www.ligo.caltech.edu/advLIGO/, http://wwwcascina.virgo.infn.it/advirgo/ 
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( |Portegies Zwart & M cMilla npOOOl |Gurkan et al.|2004t |Portegies Zwart et al.|2004l preitag et al. 
2006a] ). 



In particular, Freitag et al. ( 2006b|a ) described in detail the requirements from the point of 
view of the host cluster to form an IMBH in the center of such a system. The cluster cannot 
have "too hard" binaries, the time to reach core-collapse must be shorter than 3 Myr and the 
environmental velocity dispersion cannot be much larger than ~ 500 yr -1 . Under these conditions 
the authors find that the mass of the VMS formed is > 100M o . 

The later evolution of the VMS is not well understood, nor are the necessary conditions that 
it not evolve into a super-massive star (SMS) (see for instance |Amaro-Seoane & Spurzemp OOl 



Ama ro-Seoane et al.|2002 , Amaro-Seoane|2004] and the references in their work) in this particular 
scenario, nor are the factors that could limit the mass of such an object such that it not collapse 
into an IMBH. The process depends on a number of factors and assumptions, such as e.g. the role 
of metallicity, winds (see e.g. Belkus et al. 2007 [ though it is rather unclear how to extrapolate 
the results they obtain, which are limited to stars with masses of maximum 150M o to the masses 
found in the runaway scenario, which are typically at least one order of magnitude larger) and the 



collisions on to the runaway star from a certain mass upwards. On the other hand, Suzuki et al 



(2007) investigated the growth of a runaway particle by coupling direct N— body simulations with 
smooth particle hydrodynamics, (SPH) analyzed the evolution of the star and found that stellar 
winds would not inhibit the formation of a very massive star. More recently, Glebbeek et al 



(2009) considered the effects of stellar evolution on the runaway collision product by analyzing 
the succession of collisions from a dynamical evolution. For their low-metallicity models, the final 
remnant of the merger tree is expected to explode as a supernova, and in their high-metallicity 
models the possibility of forming an IMBH is negligible, finishing with a mass of 10-14 M Q at the 
onset of carbon burning. But as a matter of fact, these stars develop an extended envelope, so that 
the probability of further collisions is higher. Glebbeek et al. ( 2009[ ) did not change the masses 
in the dynamical simulation accordingly. In any case, self-consistent direct- summation V-body 
simulations with evolution of the runaway process are called in to investigate the final outcome. 
Henceforth we will assume that IMBHs do form; in that case, the formation of a binary of them 
(BBH) in a cluster can be theoretically explained in two different ways. 

( i) The double -cluster channel: 

In this scenario, two clusters born in a cluster of clusters such as those found in the Antennae 
galaxy are gravitationally bound and doomed to collide (see Amaro-Seoane & Freitag| 2006 for 
a detailed explanation of the process and their references). When this happens, the IMBHs sink 
down to the center of the resulting merged stellar system due to dynamical friction. They form a 
BBH whose semi-major axis continues to shrink due to slingshot ejections of stars coming from the 
stellar system. In each of the processes, a star removes a small fraction of the energy and angular 
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momentum of the BBH, which becomes harder. At later stages in the evolution of the BBH, 
GW radiation takes over efficiently and the orbit starts to circularize, though one can expect these 



systems to have a residual eccentricity when entering the LISA band ( |Amaro-Seoane & Freitag 



2006, Amaro-Seoane et al. 2010). This detector will typically be able to see systems of binaries 
of IMBHs out to a few Gpc. For this channel and volume, the authors estimated an event rate of 
4-5 yr- 1 . 

(ii) The single -cluster channel: 



Giirkan et al. (2006) added a fraction of primordial binaries to the initial configuration in the 



scenario of formation of a runaway star in a stellar cluster. In their simulations they find that not 



one, but two very massive stars form in rich clusters with a binary fraction of 10%. Fregeau et al 



(2006) investigated the possibility of emission of GWs by such a BBH and estimated that LISA 
and Advanced LIGO can detect tens of them depending on the distribution of cluster masses and 
their densities. More recently, Gair et al. (2009) recalculated the rate of Fregeau et al. ( 2006[ ) for 
the case of the proposed Einstein Telescope could see and quoted a few to a few thousand events 
of comparable-mass BBH mergers of the single-cluster channel. 



2.2. Growing up (shrinking down): The role of triaxiality on centrophilic orbits 



In the case of the double-cluster channel, the cluster, which is in rotation, results from the 
merger of the two initial clusters and may develop a triaxiality sufficient to produce sufficient cen- 
trophilic orbits. These "boxy" orbits, as seen by jBerczik et al. ( 2006[ ), are typical of systems that 
do not possess symmetry around any of their axes. In contrast to loop orbits, a characteristic of 
spherically symmetric or axisymmetric systems, "boxy" orbits bring stars arbitrarily close to the 
center of the system, since it oscillates independently along the three different axes. Therefore, 
such stars, due to the fact of being potential sling-shots, can feed the process of shrinkage of the 
BBH semi-major axes by removing energy and angular momentum from it after a strong interac- 



tion. In the strong triaxial systems of |Berczik et al.| ( [2006[ ), the rotation caused in the process of 
merger creates an unstable structure in the form of a bar. Within the bar the angular momentum 
will not be conserved and thus the BBH loss-cone is full due to stars on centrophilic orbits, inde- 
pendently of the number of stars Af*. On the other hand, for BBH systems in clusters, the role of 
Brownian motion has an impact in the loss-cone (see Amaro-Seoane et al. 2010[ ). In the models 



of Amaro-Seoane & Freitag (2006), the initial conditions are a realistic parabolic merger of two 



stellar clusters. The resulting merged cluster does not show the strong axisymmetry of Berczik 
et al.| ( |2006| ). In the simulations we address for the results of this work, the BBH (of IMBHs) is not 
stalling, in spite of the reduced number of centrophilic orbits due to the architecture of the stellar 
system. 
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Fig. 1. — Triaxiality of the resulting merged cluster for different mass fractions (upper panel) 
and the mass fraction 0.5. We calculate the semi-major axes of the ellipsoid of inertia a, b and c 
(where a > b > c) according to four different mass fractions which, in turn, are distributed on the 
basis of the amount of gravitational energy. The shorter the distance to the center of the resulting 
cluster, the lower the mass fraction. Displayed are b/a (solid lines) and c/a (dashed lines). The 
lower panel shows the shape indicators for the mass fraction 0.5, together with the evolution of the 
parameter r, an indicator for the triaxiality of the system, which tends to one as time elapses; i.e. 
the system tends to be oblate. The evolution of r is similar for the rest of mass fractions 
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In Figure [TJ we show the role of the cluster symmetry explicitly by depicting the evolution of 
the triaxiality of the cluster formed as a result of the merger of the two clusters for our fiducial 



model in the case of the double-cluster channel (which is the reference model of |Amaro-Seoane 
& Freitag 2006). After a merger which is the result of a parabolic orbit, the final system is oblate 
rather than prolate; i.e. a ~ > c, where a, b and c are the cluster axes. At the outskirts the resulting 



merged cluster is flatter and at the center the binary of IMBHs makes it rather spherical. Amaro 



Seoane et al. (2009 1 addressed the single-cluster channel scenario after the formation of the IMBHs 



and used additional simulations to further evolve the BBH. They used scattering experiments of 
three bodies including relativistic precession to 1st post-Newtonian order, as well as radiation 
reaction caused by GW, so that they did not have to integrate every single star in the cluster to 
understand the posterior evolution of the BBH. In their work, between the strong encounters, a and 



e of the BBH were evolved by resorting to the quadrupolar formulae of Peters ( 1964) 



The BBH will have completely circularized when it reaches the frequencies probed by Ad- 
vanced LIGO and the ET, because the emission of GWs takes over the dynamics of the system. 



2.3. Death 

While the emission of GWs is present all the time from the very first moment in which the 
BBH is formed, the amplitude and frequency of the waves is initially so low that no present or 
planned detector would be able to register any information from the system. Only when the semi- 
major axis shrinks sufficiently does the frequency increase enough to "enter" the LISA band, which 
we assume starts at 10~ 4 Hz. The BBH then crosses the entire detector window during its inspiral 
phase, as we can see in Figure|2j We depict the evolution of a BBH of mass 439.2 + 439. 2M . The 
reason for this particular choice of masses is to give the reader a point of reference to understand 
the whole picture. Recently |Amaro-Seoane et aT](|2010|) included the effect of the rotation of the 



host cluster and addressed the dynamical evolution of the global system with a BBH of that mass. 
The authors have shown that LISA will see the system of Figure [2] with a median SNR of few tens. 
The fact that the system merges outside its band prevents LISA from observing the loudest part 
of the BBH coalescence. In order to follow the system at this early stage of its evolution in the 
LISA band, a simple post-Newtonian approach suffices for modeling the GW radiation. We are 
far enough from the highly relativistic regime and only the inspiral phase of the BBH coalescence 
is visible to the space antenna, with a rather low SNR compared to the posterior phases in the 
evolution of the BBH. 

As the binary system depicted in Figure [2] leaves the LISA band and enters the strong field 
regime, higher order post-Newtonian corrections and eventually input from numerical relativity 
simulations need to be considered in order to model the GW waveform. Three reference frequen- 
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f(Hz) 

Fig. 2. — Amplitude of the GW emitted by a system of two equal-mass IMBHs of total mass 
878. 4M Q at 100 Mpc as seen by different GW observatories. Note that we have multiplied \h(f)\ 
by a factor 2 \ff, with / the frequency of the system. This is required in order to be able to compare 
it with the sensitivity curve of the different detectors (see section [3] for more details). From left 
to right we depict the sensitivity windows of the future space-borne LISA (dashed, gray curve), 
the Einstein Telescope (dotted, gray curve) and Advanced LIGO (solid, gray line starting sharply 
at 10 Hz). The strain of the BBH of IMBHs spends most of its inspiral in the LISA band, whilst 
the ringdown and merger occur at higher frequencies, only observable by ground-based detectors. 
Notably, the ET captures an important extent of the inspiral as well as the whole ringdown and 
merger. The averaged SNR produced by this system would be SNR LIS a = 854, SNR ET = 7044 and 
SNRacIvligo = 450. The BBH system spends approximately 0.2 yrs to go from / = 0.01 Hz (well 
into the LISA band) up to the lower cut-off frequency of Advanced LIGO, 10 Hz. These two points 
are pinpointed on the plot 
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cies in the evolution of a compact BBH that approaches its merger are the innermost stable circular 
orbit (/isco) of a test particle orbiting a Schwarzschild black hole, the light-ring frequency (/lr) 
corresponding to the smallest unstable orbit of a photon orbiting a Kerr black hole and the fun- 



damental ringdown frequency (/frd) of the decay of the quasi-normal modes computed by |Berti 
etaLl ( |2005T ). 



For the binary system shown in Figure[2j the values of these three frequencies are /isco I 878.4m© - 
5 Hz, /lr|878.4m — 14.2 Hz and /frd|878.4m q — 21.4 Hz. Should such a binary exist at a distance 
of 100 Mpc (the choice for this number is based solely on the fact that is is easily scalable), and if 
it were to be detected with Advanced LIGO, it should produce a sky-averaged SNR of ~ 450, as- 
suming a low frequency cut-off of 10 Hz. To that total SNR, the contribution of parts of the inspiral 
happening before the system reaches the characteristic frequencies /isco> Ilr an d /frd would be 
0%, 37% and 95% respectively. Figure [3] illustrates the same percentages for binaries with total 
masses between 100 and 2000 M Q . It is immediately noticed that for the binaries of IMBHs of 
interest to this study, most of the SNR that these binaries will produce in Advanced LIGO comes 
from the last stages of the BBH coalescence. 

We can estimate the time that the binary system takes to evolve from / = 0.01 Hz, a frequency 
where the BBH can be seen by LISA, to the lower cut-off frequency of 10 Hz of Advanced LIGO 
or of 1 Hz of the ET. A lower order approximation based on the Newtonian quadrupole formula 
( |Peters| [T964 ) leads to the following expression for the evolution of the frequency in terms of the 



chirp mass M. = (m l m2) 3 ^ 5 M to ^ 5 and frequency of the system 

d l = 9 l 7l s/i M yi f n/^ (1) 
dt 5 

We find a delay of only 0.2 yrs (80 days) for a BBH with total mass M = 878. 4M to go from 0.01 
Hz to the beginning of the Advanced LIGO band and almost similar numbers to the beginning of 
the ET band (the evolution of the system is extraordinarily quick in the late inspiral phase, which 
explains the fast evolution from 1 to 10 Hz). In view of these figures, LISA could be used as an 
"alarm" to prepare ground-based detectors to register the final coalescence in detail, the death of 
the BBH as such, by adjusting their sweet spots (the most sensitive part of the detector) to the 
particular BBH. The high accuracy of parameter estimation during the inspiral phase of which 
LISA is capable could be combined with the information obtained from the large-SNR triggers 
that the BBH merger and ringdown will produce in Advanced LIGO or ET and thereby achieve a 
more complete characterization of the system. 



- 10- 




Total mass (M Q ) 

Fig. 3. — Percentage of the total SNR produced by IMBH inspiralling signals cut at the three 
reference frequencies /isco;/lr an d /frd- The SNRs have been calculated using the noise curve of 
Advanced LIGO for signals placed at 100 Mpc of the detector starting at 10 Hz and with a total 
mass between 100 and 2000 M Q . Whereas the SNR computed up to /isco constitutes more than 
50% of the total SNR for systems with total mass below 100M o , it is the merger and ringdown 
parts of the coalescence (after / LR and /frd) that contribute most to the SNR as the total mass of 
the system increases above a few hundreds of solar masses 
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Fig. 4. — IMBH systems as will be seen in the time-domain output strain of the detector by the 
Advanced LIGO interferometer at the Livingston site. Three signals corresponding to equal-mass, 
non-spinning IMBH systems with total mass 400 M & < M total < 700 M & and random orientations 
and polarization angles have been placed at 1 Gpc from the detector with a starting frequency of 
10 Hz. The LI interferometer strain has been modeled by Gaussian noise colored with the design 
sensitivity curve expected for Advanced LIGO. Depending on their distance and orientation, the 
signals could be spotted by eye, which gives an intuitive idea of the kind of "bright" (in terms of 
GW emission) sources they are 
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3. BBH waveform model and expected SNR 



Accurate theoretical modeling of the gravitational radiation h(t) emitted by a BBH is key 
to improving BBH detectability and parameter estimation. While post-Newtonian (PN) theory is 
valid to model the early inspiral phase of the BBH evolution in the LISA band, an exact description 
of the merger and ringdown stages is only possible via numerical relativity (NR) calculations. 
Several approaches have been proposed to match PN and NR in order to obtain a full waveform 



for the entire coalescence, most notably for non-spinning configurations in Buonanno & Damour 
( fT999| [2000] ) ; |Buonanno etaE] ( [2009] ) ; | Ajith et al.| ( [2007] [2008] ) and recently introducing spins for 
non-precessing configurations in |Pan et a l. ( 2010}; [Ajith et al.| ( |2009[ ). 



For the purpose of the SNR and horizon distance calculations shown in this paper, we have 
chosen a new procedure for the construction of hybrid PN-NR waveforms in the frequency domain 
proposed by Santam aria et al.| ( |2010 ). The construction procedure developed for matching PN and 
NR data is sketched in Figure [5J further details on the fitting procedure can be found in Santamaria 



et al. (2010). In essence, the model consists of a phenomenological fit to hybrid PN+NR waveforms 
for spinning, non-precessing BBH systems with comparable masses (0.15 < 77 < 0.25). This is 
compatible with the simulations of Giirkan et al. ( 2006| ), which find pairs of VMSs with mass ratios 
close to 1. The phenomenological waveforms are parameterized by three physical parameters: the 
total mass of the binary M, the symmetric mass ratio r\ and the spin parameter \ = (1 + ^)Xi/2 + 
(1 — 8)xi/2, where 8 = (mi — m%)jM, Xi = Si/mj and Si is the spin angular momentum of the 
z'th black hole. Only the dominant mode £ = 2, m = 2 of the gravitational radiation enters in the 
model; the effects of the higher modes, which become increasingly significant as r) decreases, are 
neglected. For comparable-mass scenarios, this restriction does not substantially affect our results. 
The waveform model used in our calculations is valid for binary systems with aligned spins only, 
which represents a first step towards the incorporation of the spins of the BHs in full inspiral- 
merger-ringdown (IMR) models. A discussion including precessing systems with arbitrary spins to 
this date still requires work regarding availability of NR simulations and development appropriate 
techniques to match precessing merger data with PN. Work in this regard is ongoing. 

For the results presented in this article we choose to focus on three configurations: (i) equal- 
mass, non-spinning, (ii) equal-mass, equal-value spins aligned with the direction of the total angu- 
lar momentum and magnitude x = 0.75 and (iii) mass ratio 1 : 3 (jj = 0. 19), non-spinning. Provided 
with the full IMR waveforms given by our model, we are interested in assessing the detectability 
of these systems by current and future GW observatories. The SNR of a model waveform with re- 
spect to the output stream of the detector is the quantity typically quoted to signify the detectability 
of a signal. The SNR p produced by a GW signal h(t) in a detector can be computed as (see e.g. 
Thorne| 1987j |Finn|1992| ) 
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Fig. 5. — Waveform model employed in the SNR calculations of section [3] for the equal-mass, 
non-spinning BBH scenario. The amplitude of the numerically simulated and J^ + -extrapolated 
£ = 2,m = 2 mode (Reisswig et al. 2009| ) is attached to a PN calculation based on the stationary 
phase approximation that incorporates terms up to 3rd PN order. The amplitudes are stitched at a 
frequency Mf ~ 0.02 to produce a full IMR waveform. Note that the magnitudes displayed in the 
plot are dimensionless and can be scaled to account for different BBH masses 
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p=4/ -^ur df= l -^u^ dlnf > (2) 







where h(f) is the Fourier transform of the strain h(t). In the last equation S n (f) represents the 
one-sided noise spectral density of the detector. 

Figure [6] provides a graphical representation of the detectability of several BBHs by current 
and future generations of ground-based GW detectors. Displayed are the design sensitivity of 
current initial LIGO and Virgo (sensitivities that have been met or approximately met during the 
S5/VSR1 data taking), the proposed noise curves of Advanced Virgo and two possible configu- 
rations of Advanced LIGO (broadband or "base" and optimized for 30 — 3OM BBHs) and the 
designed noise budget for the Einstein Telescope in its broadband, "base" and "xylophone" (Hild 
et al.|2010] ) configurations. The hybrid waveforms constructed according to the model of Fig.[5]for 



the three chosen configurations has been conveniently scaled to represent BBH systems with total 
mass 200,500 and 1OOOM . 

As the right-hand-side of Equation[2]suggests, plotting the quantity 2 \h(f) \ yff versus a/ Sh (/) 
allows for direct visual comparison of the importance of each of the stages of the BBH coalescence. 
The three frequencies /isco, /lr and the Lorentzian ringdown frequency / LRD = 1 .2/ F rd are marked 
on top of the configuration (i) waveforms with solid circles, squares and open squares respectively. 
One can immediately appreciate that systems with total mass above 500 M fall almost completely 
below the 40 Hz "seismic wall" of the initial LIGO detectors; however they will become very 
interesting sources for the second generation of GW interferometers and the proposed Einstein 
Telescope. Indeed, as we show in Section [2| they will also be seen by the future space-borne 
LISA. Additionally it is easy to appreciate why the hang-up configuration (iii) will produce larger 
SNRs and therefore will be seen to further distances, for it merges at higher frequencies with re- 
spect to a non-spinning configuration with the same 77. In contrast, for equal spin values, the SNR 
decreases with smaller symmetric mass ratios (compare configurations (i) and (iii)). 

In Figure [7] we compute the SNR expected for these sources in each of the above-mentioned 
detectors as a function of the redshifted total mass of the system M z = (1 +z)M BBH , for optimally- 
oriented and -located sources at a distance of 6.68 Gpc (z=l). We have cross-checked our SNR 
results with those computed by |Boyle et al.| ( |2009] ) for initial LIGO at 100 Mpc, obtaining SNR 



values within 1% of those quoted by them. A direct comparison with Gair et al. (2009) shows a 



disagreement of ~ 30% in the computed SNRs for Advanced LIGO and the ET, which might be 
attributed to the different waveforms and cosmological model used. 

Unsurprisingly, the SNRs calculated for the third generation of ground-based detectors sur- 
pass the expectations for initial and Advanced Virgo and LIGO at all masses. At z = 1, SNRs of 
the order of 0.1 are expected for current LIGO and Virgo interferometers for binaries with total 
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Fig. 6. — Hybrid waveform for three BBH configurations scaled to various IMBH masses. From 
top to bottom we show BBH systems with total mass 1000,500 and 2OOM , in blue, green and 
red respectively. Solid lines correspond to the equal-mass, non-spinning configuration (i), dashed 
lines to the equal-mass, x — 0-75 configuration (ii) and dotted lines to the non-spinning, q = 3 
configuration (iii). The sources are optimally oriented and placed at 100 Mpc of the detectors. 
The symbols on top of configuration (i) mark various stages of the BBH evolution: solid circles 
represent the ISCO frequency, squares the light ring frequency and open squares the Lorentzian 
ringdown frequency (corresponding to 1 .2 times the fundamental ringdown frequency /frdX when 
the BBH system has merged and the final BH is ringing down. Currently operating and planned 
ground-based detectors are drawn as well: plotted are the sensitivity curves of initial LIGO and 
Virgo, two possible configurations for Advanced LIGO (zero detuning and 30 — 30M o BBH opti- 
mized), Advanced Virgo and the proposed Einstein telescope in both its broadband and xylophone 
configurations 
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Fig. 7. — Signal-to-noise ratio as a function of the redshifted total mass of the BBH for the present 
and future generations of GW detectors and LISA. The sources are placed at a distance of 6.68 Gpc 
(z = 1) and the SNRs correspond to sources optimally oriented and located. Solid lines indicate 
SNRs for the equal-mass, non-spinning configuration (i); for Advanced LIGO and ET we have 
included the SNRs produced by configurations (ii) and (iii) as well, indicated with dashed and 
dotted lines respectively 
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mass up to a few hundred M & . These values would scale to SNRs of ~ 1 at the closer distances 
typically surveyed by the initial interferometers, which are of the order of tens to few hundreds of 
Mpc. Signals with single-interferometer SNRs below the commonly-used threshold of SNR=5.5 
are most likely to be missed, therefore IMBH binary coalescences are of limited interest for the 
first-generation detectors -and even for their current, enhanced counterparts-, which are in turn 
most sensitive to neutron star binaries and stellar-mass back hole binaries. 

Advanced LIGO and Virgo will be able to measure averaged SNRs of the order of 1 — 20 at 
Z — 1, with a maximal response to BBH systems with total mass in the range of 300 to 800 M Q . 
Therefore, should binaries of IMBHs exist in our neighboring Universe, the second-generation of 
GW interferometers should be able to detect them, for they will be loud enough to stand above 
the detector noise. As we will see later, the non-negligible IMBH binary coalescence rate for 
advanced LIGO indicates that these sources should be taken into account in future match-filtered 
searchers for inspiral binaries. Due to the moderate SNR of the signals, a potential detection might 
be sufficient to confirm the existence of binaries of IMBHs but not enough to determine all its 
parameters, such as mass ratio and spins, with sufficient accuracy; For the advanced detectors, 
Cutler & Flanagan ( 1994) estimate an accuracy of ~ 1% in the reduced mass ji = mim 2 /M using 



PN templates with negligible spins at SNR~ 10. When they take the spins into account, A/x 
increases by ~ 50. It is expected that the use of full IMR templates will improve these figures, 
but the final accuracy will always be limited by the slightly-above-threshold SNRs expected in 
Advanced LIGO and Virgo. 

The Einstein Telescope will instead measure SNR values within the 10 2 range at z ~ 1, and it 
is expected to be sensitive to binaries with total masses of the order of 1O 4 M , a significantly larger 
range than that surveyed by Advanced LIGO and Virgo. It is noticeable how the ET xylophone 
configuration increases the detectability of binaries with masses above 1000 M Q with respect to the 
broadband ET configuration. This is due to its improved sensitivity precisely at frequencies in the 
range of 1 — 30 Hz, which is where systems of mass above thousands of solar masses accumulate 
most of their SNR (see Figure [6]). The large-SNR events that binaries of IMBHs would produce 
in the ET are of great importance for astrophysics; the reason being that as the SNR increases, the 
accuracy in parameter estimation also does. In the limit of large SNR, the inverse of the Fisher 
information matrix T a b = (h. a \hj,) is an estimator of the errors in the recovered parameters. At the 
relative large SNRs potentially produced by IMBH binaries in the ET, the possibility of extracting 
their mass ratios and spins to high accuracies would be revolutionary for characterizing the IMBHs 
populations. 

As for LISA, BBHs with masses of hundreds of solar masses will be seen with a moderate 
SNR (see Amaro-Seoane et al.|2010 ,for a detailed study of the parameter extraction) - it is only at 



masses above tens of thousands of solar masses that LISA will start taking over the ground-based 
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observatories. Although the space antenna will be most sensitive to BBH binaries with masses 
in the range of 1O 6 ~ 7 M , the possibility that it can act as a complementary observatory for the 
Einstein Telescope for IMBH binaries is very promising. Parameter accuracy studies for IMBHs in 
LISA are already available using the inspiral part of the coalescence (including also non-negligible 
eccentricities; see Amaro-Seoane et aL]|2010[ ), and indicate that masses and sky positions will be 



recovered with a high accuracy level. In order to complete the characterization of IMBHs with the 
information given by the second and third generations of ground-based detectors, a comprehensive 
study of parameter recovery taking the BBH coalescence into account is very much desirable. 



4. Event Rates 



Miller (2002) estimated for the first time the event rate for intermediate-mass mergers of 
IMBHs (i.e. typically stellar black holes merging with IMBHs) in clusters by calculating the 
luminosity distance for the inspiral, merger and ringdown (Flana gan & Hughes|19 98) out to which 
these three stages can be detected with an SNR larger than 10. In his approach, the maximum 
distance for the detector was 3 Gpc (z ~ 0.53), with no cosmological corrections. The event rate 
was calculated as R = J (4/3) D(M) 3 u{M) n ng f(M) dM. In this equation n ng is the number density 
of globular clusters, which was taken to be n n „ ~ 8/z 3 /Mpc 3 , as in the work of Portegies Zwart & 



McMillan (2000). The rate of coalescence of stellar-mass compact objects with the IMBH is v{M) 



and f(M) = dN /dM is the mass distribution of massive enough black holes in clusters. Obviously, 
/ f{M) dM = / tot < 1. |Miller|([2002|) uses the estimation of|Flanagan & Hughes] (|1998|) for it and 



finds that a few per year should be detectable during the last phase of their inspiral. Two years later, 



Will (2004) revisited the problem using matched filtering for the parameter estimation, an updated 
curve for the sensitivity of the detector and restricted post-Newtonian waveforms to calculate an 
analytical expression for the luminosity distance D L ; his estimation is a detection rate for binaries 
of about 1 per Myr in a mass-range of (10 : 100)M o . 



Subsequently, Fregeau et al. (2006) calculated the number of events that initial and Advanced 
LIGO (and LISA) could see from the single-cluster channel. In their estimation, they assume 
that the VMSs formed in the runaway scenario do not merge into one, but evolve separately and 
eventually each form an individual IMBH, following the numerical results of the Monte Carlo 



experiments of Giirkan et al. (2006). They derive a generalized form of the event rate which can 
be summarized as follows 



R 



dN„ 



dtn 



d 2 M. 



SF 



dV c dt e ' 



dt e dV c 
dto dz 



dNl 



dM SF . cl dM cl 



dM c{ dz. 



(3) 
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In this expression, d 2 M SF /dV c dt e is the star formation rate (SFR) per unit of comoving volume per 
unit of local time; g cl is the fraction of mass that goes into the massive clusters of interest; g is the 
fraction of massive clusters which form IMBHs; dt e /dt = (1 +z) _1 is the relation between local 
and observed time; dV c /dz is the change of comoving volume with redshift; and dN^/dM SFc idM c i 
is the distribution function of clusters over individual cluster mass M c \ and total star-forming mass 
in clusters M SFc i; z max is the maximum redshift that the detector is capable of observing. 

In order to compute the integral above, an estimation of the observable volume of each indi- 
vidual detector is required. A commonly-used measure of the reach of a detector is the horizon 
distance Dh, defined as the distance at which a detector measures an SNR = 8 for an optimally- 
oriented and optimally-located binary, i.e. an overhead, face-on orbit. Suboptimally located 
and oriented sources are detected with SNR = 8 at closer distances. However, in order to com- 
pute merger rates, an average distance over all possible orientations is more meaningful. The 
orientation- averaged distance represents the cube root of the total volume to which a detector is 
sensitive, assuming uniformly-distributed sources, and is 2.26 times smaller than the horizon dis- 
tance D h ( |Finn & Chernoff||1993| ). Moreover, at the distances that Advanced LIGO and the ET 



are expected to survey, a certain cosmological model needs to be assumed. We adopt the standard 
ACDM universe with parameters given by the first five years of the WMAP sky survey (Hin 



shaw et al.||2009| ). These are fi A = 0.73, Vl h = 0.046, Q c = 0.23, H = 70.5kms _1 Mpc _I and 
to = 13.72 Gyr. Using the full IMR waveforms described in Sec.|3]and the corresponding redshift 
function z(d) for the ACDM model, we compute the orientation-averaged distance D/,/2.26 for 
non-spinning systems with symmetric mass ratio rj = 0.25, 0. 19 and for an equal-mass system with 
spins aligned in the direction of the angular moment and total spin x — 0.75 (hang-up configu- 
ration). Computing Dh(M z ) is essentially equivalent to inverting Eq. [2] for p(d,M(d)) for p = 8 
using the relation z(d) given by the cosmological model. The results for D h {M z ) for our three 
configurations can be seen in Fig. [8] for Advanced LIGO and ET. 

The limit for the integral in z is thus given by the observable volume of the detector of interest, 
quantified using the orientation- averaged distance. The maximum redshift z max is in turn a function 
of the mass of the binary system, its configuration -mass ratio, spin- and the particular waveform 
model used in the calculation of the horizon distance. We present our values for z max in Fig. |9} 
which essentially shows the same data as Fig.[8]but expressed in terms of redshift versus total mass 
of the binary. As we see in Fig. [8j the maximum values for the orientation-averaged distance for 
ET are as large as z ~ 10 for configuration (ii) at M z ~ 3OOM . This implies that the ET will be 
able to probe the different proposed scenarios to produce the first generation of black hole seeds, 
as pointed out by |Sesana et al.| ( [2009] ). However, at these large cosmological distances, the stellar 



formation rate is unknown and the validity of the rate integral cannot be stated. We therefore set a 
maximum cutoff value of z max = 5 in the computation of Eq.[5j The final value of Zmax C^bbh) that 
we have used in the computation of the rates is show in figure [9] for our three particular physical 
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Fig. 8. — Orientation-averaged distance versus redshifted mass for three binary configurations 
obtained with the design sensitivity curves of Advanced LIGO and the Einstein Telescope. The 
solid, dashed and dotted lines corresponds to the configurations denoted in the text as (i), (ii) 
and (iii) respectively. Note the ~ 40% increase in reach given by the hang-up configuration with 
X = 0.75 with respect to the non-spinning case 
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configurations. 

Regarding the term in the dM ci integral in [5J three different parameterizations of the SFR are 
available in the literature, i.e. /?sfi,2,3 as summarized by equations 4, 5 and 6 of IPorciani & Madau 



(2001 ). The three models are similar for distances up to z ~ 2, where SFR peaks, differing from 
there on. For the results shown here, we have compared the three of them. As for the distribution 
of cluster masses, the factor can be approximated as 



dNl 



f{M c 



dM SF , ci dM cl f M c i /(M c i) dM cl ' 



(4) 



In order to compute the integral in the denominator, we take dN c \/dM c \ oc 1/M^, following the 
power law form observed by Zhang & Fall ( 1999[ ) for young star clusters in the Antennae. The 
validity of assuming the same law for the larger volume of the Universe surveyed by Advanced 
LIGO or the ET is, however, a generalization not based on direct observations. Thus, we should 
take this premise with care. By assuming an efficiency factor of /gc ~ 2 x 10~ 3 , based in the 
results of Giirkan et al. (2006), we can set the values for the limits in M cl = /gc/^bbh according 



to the masses of the IMBH binaries of interest. In our case, taking the standard definition of IMBH 
into account, this is M chmax /M chmin = M bbh. max /M bbh, min = 2 x 10 4 M o /2 x 1O 2 M . The integral 
can now be expressed as 
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GC 



hi(M cl , max /M cl 



Mbbh., 



dM ] 
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Mbbh., 



-max (Mbbh) 
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1 dV r 



dz 



dz, (5) 



where A/BBKmaxOin) is the range of total mass of the BBH that we are considering, SFR t (z), i 
1, 2, 3 is any of the three considered stellar formation rates of |Porciani & Madau] ( |2001[ ) and 



Hz) 



(l+ Z )3/2 



(6) 



is the factor that relates the stellar formation rate function in different cosmologies with respect to 
the Einstein-de Sitter Universe. 

Regarding the values of the terms g and g c i, it is unfortunately very little what we know 
about the initial cluster conditions required to form an IMBH binary. These factors have therefore 
large uncertainties. Following the existing literature, we leave the fraction of massive clusters that 
form IMBHs, g, as a parameter and set it to 0.1 as an example. Nevertheless, as proven in the 



simulations of Freitag et al. (2006b I, it could be as large as 0.5. As for the fraction of mass going 
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into massive clusters, g ch previous works have also tentatively set it to 0. 1 . Nevertheless, the results 
of McLaughlin (1999) seem to indicate that this value might be too optimistic. McLaughlin (1999) 



estimates an empirical cluster formation efficiency by mass of e = M^JMq^, where M^ s is the 
initial globular cluster population and M^ s is the initial reservoir available in the protogalaxy. He 
then infers a universal value of e ~ 0.25% after evaluating a sample of 97 giant ellipticals, brightest 
cluster galaxies and faint dwarfs. He finds identical results for the Population II spheroid of the 
Milky Way. Therefore we have chosen to set our g ci to 1/400, noting that, in any case, the final 
rates scale trivially with both g and g d . 

We have evaluated the integral [5] for the three stellar formation rates we consider and the three 
binary configurations (i), (ii) and (iii) for Advanced LIGO and the ET. Our results are summarized 
in table [T] We find that the event rate does not depend strongly on the assumed stellar formation 
rate; the differences are of the order of ~ 5% for Advanced LIGO and ~ 30% for the ET. In the 
case of Advanced LIGO, the insensitivity to the stellar formation model is easily explained by the 



fact that the three models of Porciani & Madau (2001 1 are very similar until z ~ 2; our integration 



in z stops at z ~ 1 (see red curves in Fig. [9]), therefore ignoring contributions at higher z. For the 
ET, the differences among SFRi (z) , i = 1 , 2, 3 between z ~ 2 and 5 are very much attenuated by the 
rapid decrease in comoving volume at high redshift. In order to assess the effect of stopping the 
integration at z = 5, we have also computed the rates for the ET without imposing this condition. 
The differences amount to ~ 10%, therefore we conclude that the rates estimations are not very 
sensitive to very high-z SFRs. Nonetheless, we should stress again that the lack of stellar formation 
data at these large cosmological distances makes assumptions at z > 5 largely speculative and thus 
we restrict to presenting the rates calculated with z max clamped at a value of 5. Given the inherent 
uncertainties present in our approach, which can amount to a few orders of magnitude, the factors 
associated to the choice of SFR do not represent a major source of error. We therefore quote the 
results found for SFR 2 (z) only, corresponding to the star formation rate that increases up to z ~ 2 
and keeps constant afterwards. 

The event rates do, however, strongly depend on the spins and mass ratio of the binary. As 
expected, "loud" configurations like the hang-up case increase the event rate by a factor of ~ 3 
in the case of Advanced LIGO. This effect is more notable the larger the total spin of the binary 
is; the correct determination of the spin distribution of IMBH binary systems would be extremely 
useful to further quantify the impact on the total rates. Smaller mass ratios decrease the event 
rates, the difference between 77 = 0.25 and i] = 0.19 being also approximately of a factor of 3. 
The differences are not so extreme in the case of the ET, due to the fact that we are cutting off 
z max at a value of 5 and, thus, neglecting contributions at higher redshift which might amount to 
another ~ 10%, as discussed above. Even so, it is evident that the expected rates increase with 
the total spin \ an d the symmetric mass ratio i]. The exact computation of the total rates for all 
possible physical configurations of the binary would necessitate further integration on the mass 
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Fig. 9. — Orientation-averaged distance expressed in terms of maximum redshift z max for the 
Advanced LIGO and the ET detectors up to z = 5 versus intrinsic total mass of the BBH. Red lines 
are for Advanced LIGO and blue are for the ET. The solid, dotted, dashed curves represent the 
same physical configurations displayed in figure [8] 
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Table 1: Event rates formed in the single-cluster channel for IMBH binaries potentially observable 
by Advanced LIGO and the ET per year. We take g = 0.1, g c i = 1/400 as standard scaling values 



Detector Configuration R ( — ) — ; 

5 [\0.lJ \ 1/400 

Advanced LIGO r? = 0.25, % = °- 85 

rj = 0.25, x = 0-75 2.36 

77 = 0.19, x = 0-31 

Einstein Telescope 77 = 0.25, x = 21.5 

rj = 0.25, x = 0-75 24.0 

77 = 0.19, x = 16.8 



ratios and spins of the system. At present it is not clear what that distribution might be, therefore 
we simply summarize our results for all configurations under consideration in table [T] We find 
event rates of the order of ~ 1 (g/0.1) (g c i/l/400)yr _1 for Advanced LIGO and of the order of 
~ 20(g/0.1) (g c i/l/400)yr _1 for the ET. These rates assume formation of IMBH binaries in the 
single-cluster channel. 



For Advanced LIGO and a non- spinning configuration, Fregeau et al. (2006 ) found a expected 
detection rate of 10(g/0.1) (g c i/0.1)yr _1 (0.25 yr -1 in our units), adopting a reach of 2 Gpc (z max ~ 
0.37). This estimation, calculated assuming a ringdown-only search, underestimates the volume 
of the Universe that Advanced LIGO would observe using full IMR templates, as we show in 
Fig. [8] Our results for non-spinning systems are compatible with this observation; we quote rates 
larger than those of Fregeau et al. ( 2006| ) by a factor of 3 but otherwise of a compatible order of 
magnitude. In the case of the ET,|Gair et al](|2009|) quote a rate of ~ 500 (g/0.1) (g c i/0.1)yr -1 (~ 



13yr _1 in our units) for non- spinning configurations, a value of the same order of magnitude than 
the one found by us. The fact that Gair et alT] ( |2009] ) use a different family of waveforms to model 
the coalescence and a fitted formula for the the averaged distance based on EOBNR waveforms, 
together with the slightly different values in the integration limits explain the discrepancies in 
the exact figures. Therefore, our new results for the rates of IMBH binary coalescence for non- 
spinning systems agree reasonably well with previous works that used similar methods but different 
waveform and detector models. In addition, we have now quantified the effect that the spins of the 
BHs will have on the expected rates and have found it to be of a factor ~ 3 for total spins of the 
binary as high as x = 0.75 



Thus far we have concentrated on the single-cluster channel scenario. 



Amaro-Seoane & Fre- 
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itag (2006) give a prescription to calculate an estimate of the event rates for the double-cluster 
channel. This was based in the fact that the only difference, in terms of the event calculation, be- 
tween both astrophysical scenarios involves, firstly, that in the double-cluster channel there is one 
single IMBH in one cluster and, secondly, that these two clusters must collide so that the IMBHs 
form a BBH when they sink to the center due to dynamical friction. As explained in section 4 



of Amaro-Seoane & Freitag (2006), the connection between the event rate estimation of the two 
channels is 

r doub = p merg gr 8 ^, (7) 

where r doub is the event rate of the double-cluster channel, T sing of the single-channel and P merg is 



the probability for two clusters to collide in the scenario of Amaro-Seoane & Freitag (2006). They 
find that the unknown parameter P me rg could have values in the range P merg G [0.1, 1]. The total 
event rate assuming that both formation channels are possible would be 



i sing 



gel 



0.1/ V 1/400/ V merg V0.1 



hrr , (8) 



Assigning parameter P merg its pessimistic and optimistic limits, we can compute the lower and 
upper limit of the event rates for Advanced LIGO and the ET, assuming contributions of the two 
channels. Taking into account the double-cluster channel increases the rates upper limit by a factor 
of 2. For instance, for the equal-mass, non-spinning case, the values are 

rSxroo e [(0)0.94, 1.7] (^) (^) yr- (9) 

e [(0) 23.7, 43.0] (A.) (jj^j yr- 1 (10) 



Note that the results are quoted in terms of g and g c \; the already-mentioned uncertainty in their 
values could increase these rates by at least one order of magnitude. These event rates are encour- 
aging to address the problem of detection and characterization of systems of IMBH binaries with 
ground-based GW observatories. On the other hand, one should bear in mind that the existence of 
IMBHs altogether has not yet been corroborated, so that the pessimistic estimate is still somewhat 
"optimistic". This is why we have added a (0) in the previous rates as the absolute lower limit. 



5. Conclusions 

Even though we do not have any evidence of IMBHs so far, a number of theoretical works have 
addressed their formation in dense stellar clusters. If we were to follow the same techniques that 
have led us to discover the SMBH in our own Galaxy, we would need the Very Large Telescope 
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interferometer and next-generation instruments, such as the VSI or GRAVITY, which should be 
operative in the next ~ 10 yrs. An alternative or even complementary way of discovering IMBHs 
is via their emission of GWs when they are in a BBH system. 

The identification and characterization of these systems rely on accurate waveform modeling 
of their GW emission, which has been made possible by the success of numerical relativity in 
simulating the last orbits of the BBH coalescence and the coupling of these results to analytical 
post-Newtonian calculations of the inspiral phase. We use a PN-NR hybrid waveform model of 



the BBH coalescence based on a construction procedure in the frequency domain (see Santamaria 
et al.|2010[ for details). 



Using this hybrid waveform, we have estimated the SNR corresponding to the current and 
Advanced LIGO and Virgo detectors, the proposed ET and the space-based LISA at a distance of 
z = 1, i.e. 6.68 Gpc. The results indicate that IMBH binaries will produce SNRs sufficient for 
detection in advanced LIGO and Virgo and notably larger SNRs in the ET, thus making them inter- 
esting sources to follow up on. Eventual observations of IMBH binaries with future ground-based 
detectors could be complementary to those of LISA, which is expected to detect these systems with 
moderate SNRs and to be more sensitive to SMBH binaries instead. More remarkably, in principle, 
if LISA and the ET are operative at the same time, they could complement each other and be used 
to track a particular event. 

Furthermore, we have revisited the event rate of BBHs for various detectors and find encour- 
aging results, within the inherent uncertainties of the approach. Our estimations are consistent 
with previous works, and additionally we have quantified the effect of the total spin of the binary 
in the expected event rates. We have estimated the distance to which Advanced LIGO and the 
ET will be able to see binaries of IMBHs. This quantity depends strongly on the mass ratio and 
spins of the binary. For Advanced LIGO, equal-mass, non-spinning configurations of observed 
total mass ~ 2OO-7OOM can be seen up to z ~ 0.8. If the spins are aligned with the total an- 
gular momentum and significant (xi,2 ~ 0.75), the reach increases to z ~ 1.5 for observed total 
masses of ~ 400 M . The ET will be able to explore even more remote distances, reaching to 
Z ~ 5 and further. Our present knowledge of stellar formation at such large redshifts is incomplete, 
therefore we have computed the event rates for the ET integrating only until z = 5. We have com- 
pared three star formation models and three different configurations of the binary. The effect of 
the particular formation model is not very significant; remarkably, the physical configuration does 
however influence the final rates. We provide the results for three physical configurations, taking 
into account both single- and double-cluster channels in the binary formation. For a fully correct 
calculation of the event rate integrated over all possible configurations, more detailed knowledge 
of the distribution of spins and mass ratios of IMBH binaries formed in globular cluster is required. 

Advanced ground-based detectors are designed to be able to operate in different modes so that 
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their sensitivity can be tuned to various kinds of astrophysical objects. Considering the importance 
of an eventual detection of a BBH, the design of an optimized Advanced LIGO configuration for 
systems with M ~ 1O 2-4 M would be desirable to increase the possibility of observing such a 
system. In case an IMBH binary coalescence were detected, the recovery and study of the physical 
parameters of the system could serve to test general relativity and prove or reject alternative theo- 
ries, such as scalar-tensor type or massive graviton theories. The direct identification of an IMBH 
with GWs will be a revolutionary event not only due to the uncertainty that surrounds their exis- 
tence and their potential role in testing general relativity. The information encoded in the detection 
will provide us with a detailed description of the environment of the BBH/IMBH itself. 

We are indebted to C. Cutler, B. Sathyaprakash and B. Krishnan for discussions about the 
work. We thank C. Reisswig and D. Pollney for providing us with their Cauchy-extrapolated data 
and F. Ohme for helping us to implement the post-Newtonian formalism used for the waveform 
modeling. We are grateful to E. Robinson and D. Vanecek for comments on the manuscript. PAS 
work was partially supported by the DLR (Deutsches Zentrum fur Luft- und Raumfahrt). LS has 
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